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\ 

\i 

This  interim  report  describes 

(i)  Aerosol  particle  number  and  size  distributions  for  ambient  air  at  Mace 
i;ead  Atmospheric  Research  Station,  Carna  in  the  Nest  of  Ireland  about  50 
miles  west  of  Galway  City. 

(ii)  Intercomparison  of  particle  sizing  instrumentation  used  by  University 
College  Galway  and  'uMIST,  Manchester. 


(iii)  Calculations  of  optical  extinction  and  lidar  backscatter  to 
extinction  ratios  at  predominantly  laser  wavelengths  from  the  visible  to 


middle-IR  -  based  on  measured  aerosol  particle  size  distributions. 


(i)  Aerosol  particle  number  and  size  distribution  measurements  at  Mace 
Head. 

Introduction 

Use  was  made  of  a  Particle  Measuring  Systems  (PUS)  Active  Scattering  Probe 
-  ASASP-X  probe.  The  probe  was  realigned  and  calibrated  for  size  response 
by  means  of  polyvinyltoluene  latex  and  polystyrene  latex  monodisperse 
spherical  particles.  The  size  calibration  agreed  to  within  plus/minus  one 
channel  resolution  of  the  manufacturers  response,  and  is  given  in  Table  1 
be low . 

Table  L.  Size  calibration  of  Particle  Measuring  Systems  (PMS)  ASASP-X 
light  scattering  counter. 


Range  3  Range  2  Range  1  Range  0 

Particle  diameter 

intervals  (pm)  0.09-0.20  0.15-0.30  0.24-0.R4  0. 6-3.0 

Field  measurements  have  been  carried  out  spanning  both  1987  and  the  present 
year  at  about  one  monthly  intervals.  A  listing  of  the  field  test  dates  is 
given  in  Table  2. 


Table  2. 


Field  Measurement  Tests 


Date 

Prevaling  Wind 
Direction 

Comments 

December 

10  -  11,  1987 

Easterly 

2  s  data  throughout. 

Solenoid  valves  in  operation. 

December 

16  -  17,  1987 

Southerly 

Some  measurements  bypassing 
solenoid  valves  were  made. 

January  13  -  15,  1988 

South-westerly- 

westerly 

2s  data. 

February 

13-19,  1988 

Intercalibration  of  particle 
sizing  instrumentation  at 
UMIST,  Manchester. 

March  10 

-  12,  1988 

Westerly 

4  s  data,  no  solenoid  valves 
in  use. 

April  19 

-  22,  1988 

Southeasterly  to 
southwesterly 

Parallel  measurements  at 

Mace  Head  and  at  Great 

Dun  Fell. 

The  ambient  aerosol  field  measurements  of  December  10-11,  1937  were  made 
with  a  normally  open  solenoid  valve  in  line  to  the  ASASP-X  probe. 

Laboratory  tests,  using  ammonium  sulphate  polydispersion  of  particles  as  a 
test  aerosol,  were  carried  out  to  compare  particle  concentration  passing 
through  the  normally  open  valve  with  that  passing  directly  to  the  sizing 
probe.  Results  of  this  comparison  yielded  averaged  losses  of  12.1%,  14%, 
19.6%  and  62.8%  for  ranges  3,  2,  1,  0.  These  percentage  losses  are  mainly 
due  to  the  presence  of  a  sharp  L  shaped  bend  in  the  normally  open  valve. 
Ambient  aerosol  particle  concentration  measurements  made  on  December  10-11, 
1987  were  corrected  by  applying  the  above  laboratory  loss  correction  valves 
to  the  data.  The  rest  of  the  ambient  aerosol  data  presented  was  obtained 
via  direct  access  to  the  ASASP-X  probe. 


The  aerosol  particle  size  data  was  taken  at  2  s  intervals  for  the  first 
three  field  data  sets.  Subsequent  field  data  is  taken  at  a  4  s  interval. 
Particle  concentration  is  measured  directly  since  the  PMS  ASASP-X  is 
totally  plumbed  i.e.  each  particle  is  counted  as  it  traverses  the  laser 
beam.  Hence  particle  concentration  is  known  as  accurately  as  the  flow 


rate  through  the  particle  counter, 
were  performed  in  the  laboratory. 


Periodic  flow  calibration  measurements 
A  flow  rate  of  1  cm  s  ^  was  maintained 


for  all  ambient  aerosol  measurements. 


Analysis  of  the  Aerosol  Field  Data 

The  aerosol  particle  data  is  presented  here  in  two  formats: 

a 

(i)  Aerosol  concentration  per  era  as  a  function  of  particle  radius. 

a 

(ii)  Aerosol  concentration  per  cm  per  micron  as  a  function  of  particle 
radius. 

Format  (1)  tends  to  show  'spikes'  which  are  due  to  overlapping  size  ranges 
with  larger  number  concentration  than  in  neighbouring  size  channels  of  the 


adjacent  (lower)  size  range.  These  spikes  are  smoothed  out  using  format 
(li)  where  the  particle  concentration  is  weighted  against  the  width  of  the 
particle  channel. 

The  data  is  presented  here  in  chronological  order.  Figure  1  and  2  show 
the  uncorrected  data  of  aerosol  particle  concentration  per  cc,  per  cc  per 
micron  as  a  function  of  particle  radius  for  December  11,  1987.  The 
corrected  data,  allowing  for  particle  loss  through  the  solenoid  valve  is 
shown  in  Figures  3  and  4.  A  relatively  high  particle  concentration  is 
noted  for  the  easterly  wind  direction  -  corresponding  to  the  source  of  the 
aerosol  coming  from  the  continental  European  landmass. 

Additional  ambient  aerosol  concentration  data  is  presented  in  Fig.5  and  6 
for  Dec.  17,  1987  at  Mace  Head,  for  a  predominantly  southerly  airstream, 
Fig.  7,  Fig.  8,  Fig.  9  and  Fig.  10  shows  particle  size  data  for  five  minute 
periods  on  January  14th  and  15th,  1988. 

An  intensive  measurement  period  was  also  carried  out  for  March  lOth-llth, 
1988.  Data  is  presented  firstly  for  four  consecutive  cycles  of  duration 
twenty  five  minutes  over  the  time  period  1805  through  2120.  Figures  11, 
and  12  give  consecutive  cyclic  particle  concentration  as  a  function  of 
particle  radius  for  both  formats.  An  average  of  the  four  cycles  data  is 
given  in  Figs.  13  and  14.  Similar  data  plots  are  shown  in  Figs.  15,  16, 

17  and  13  for  five  consecutive  data  cycles  of  duration  twenty  five  minutes 
over  the  time  period  2230  through  0235,  on  March  10th  and  11th,  1983.  A 
more  comprehensive  analysis  of  these  data  sets  will  be  performed  for  the 
final  report.  In  addition,  analysis  will  also  be  made  of  the  data  set  of 
April  19-22,  1983  which  forms  a  major  part  of  the  joint  trajectory  study 
carried  out  in  parallel  with  the  UMIST  research  group  under  the  supervision 
of  Dr.  T.  Choularton,  at  Great  Dun  Fell,  over  the  same  time  period. 


(ii)  Intercomparison  of  particle  sizing  instrumentation  used  by  University 
College  Galway  and  UMIST,  Manchester. 

The  PMS  ASASP-X  probe  was  brought  to  the  Department  of  Pure  and  Applied 
Physics,  UMIST,  Manchester  on  February  13tn  1988.  Intercomparison  was 
performed  between  the  Galway  ASASP-X  probe  and  the  UMIST  ASASP  probe.  The 
inter-range  compatibility  of  the  two  probes  is  shown  in  Table  3.  The 
Galway  ASASP-X  was  adopted  as  the  standard  probe  since  it  measured  particle 
number  concentration  absolutely.  Laboratory  generated  polydispersion  of 
ammonium  sulphate  were  used  as  the  test  aerosol.  Sequential  measurements 
of  particle  number  concentration  of  the  ammonium  sulphate  aerosol  were  made 
using  the  ASASP-X  and  the  ASASP  probes.  A  detailed  analysis  of  this 
comparison  is  ongoing  and  will  be  presented  in  the  final  report. 


(iii)  Calculation  of  optical  extinction,  and  lidar  backseat  ter  based  on 
the  measured  particle  size  distributions 

A  modified  version  of  the  Mie  program  of  Dave  (1963)  was  used  for  E-0 
scattering  calculations  in  this  work,  under  the  simplifying  assumption  that 
all  particles  are  homogeneous  spheres. 


The  transmission  of  a  monochromatic  laser  beam  through  an  assemblage  ot 
particles  is  derived  from  the  Lambert-Bouguer  law.  The  extinction 


coefficient 


tfe(m  A)  is  determined  from  the  measurement  of  the  transmitted 
signal  I  after  traversal  through  the  assemblage  of  particles  of  path  length 
L: 


I 


I0exp.(-aeL) 


(1) 


TABLE  3 


The  inter-ranae 

comoati bi 1 l tv  of 

the  ASASX  and  ASASP 

ASAS-X  Ranqe  0 

=  ASASP  Ranae  0 

Ranae 

particle  size 

channnel 

parti cle  size 

ASAS-X 

ASASP 

0. 600-0. 759 

ch  1 

0. 240-0. 280 

ch  1 

0.750-0.919 

ch2 

0.280-0.320 

ch2 

0.920-1 .079 

ch3 

0. 320-0.360 

ch3 

1 . 080-1 . 239 

ch4 

0. 360-0. 400 

ch4 

1 . 240-1 . 399 

ch5 

0. 400-0. 440 

ch5 

ch  t 

1 . 400-1 . 559 

ch6 

0. 440-0. 4H0 

c  h  6 

ch2 

1 . 560-1.719 

ch7 

0. 480-0.520 

ch7 

ch3 

1 . 720-1 . B79 

ch8 

0.520-0. 560 

chB 

c  1)4 

1 . 880-2. 039 

ch9 

0. 560-0. 600 

ch9 

c  h5 

2.040-2. 199 

ch  1 0 

0.600  0.640 

ch  1  0 

c.h6 

2.200-2. 359 

ch  1  1 

0.640-0.680 

ch  1  1 

ch 

2.360-2.519 

ch  1  2 

0.680-0.720 

ch  1  2 

chB 

2.520-2.679 

ch  1  3 

0. 720-0. 760 

ch  1  3 

ch9 

2.680-2. 839 

ch  1  4 

0.  760-0.800 

chi  4 

c  h  1 0 

2.840-3.000 

ch  1 5 

0.800-0.840 

ch  1  5 

ch  1  1 

0.840-0. 880 

ch  1  2 

0. 880-0. 920 

ch  1  3 

0.920-0.960 

ch  1  4 

O. 960-1 . 000 

ch  t  5 

ASAS-X  Range  2  =  ASASP  Range _ 3 

particle  size _ channel_ _ 


0. 15-0. 16 

chi 

0.16-0.1/ 

ch2 

0. 17-0. 18 

ch3 

0. 18-0. 1 9 

ch4 

0. 19-0. 20 

ch5 

0.  20  -0.  21 

chi 

0. 21-0.22 

ch  7 

0. 22-0. 23 

chB 

0.23-0. 24 

ch9 

0.24-0. 25 

ch  1  0 

0. 25-0.26 

ch  1  1 

0. 26-0. 27 

ch  1  2 

0.27-0. 28 

ch  1  3 

0.28-0. 29 

ch  1  4 

0. 29-0. 30 

ch  15 

ASAS-X  sample  volume  1  cc/sec 

ASASP  sample  volume  0.105  cc/sec 

-  multiplication  factor  =  9.52 


where  I  is  the  initial  intensity  of  the  laser  radiation  in  the  absence  of 
particles. 

The  volume  extinction  and  backscatter  coefficients  and  0^  of  a 
polydispersion  of  spherical  particles  characterized  by  a  size  distribution 


n(r)  and  refractive  index 

m  are  given  by 

a  = 
e 

/  irr2Qen(r)dr, 

(2) 

rb  = 

1  2 

i - /  7rr^Gn(r)dr , 

4  r 

(3) 

where  Qe(ra,x)  is  the  Mie  efficiency  factor  for  extinction  for  a  particle 
with  refractive  index  m  and  size  parameter  x  =  2~nr/\  ,  and  G(ra,x)  is  the 
backscatter  gain  defined  as  4m  times  the  ratio  of  the  backscatter 
differential  cross  section  to  the  geometric  area. 


Refractive  index  values  were  chosen  to  represent  the  more  commonly  found 
atmospheric  constituents  for  maritime  aerosol:  sodium  chloride  and 
ammonium  sulphate.  Refractive  index  values  for  these  two  constituents  for 
visible  to  middle  infrared  wavelengths  are  given  in  Table  4.  Also  shown 
in  Table  4  are  calculated  values  of  total  extinction  coefficient  (m-^) 
based  on  the  averaged  measured  particle  size  distribution  data  of  March 
lOth-llth,  1988  over  the  time  period  2230-0235  (see  Fig. 18). 


Table  4.  Extinction  coefficient  (ta~1')  for  visible  up  to  middle  IR 
wavelengths  based  on  measured  particle  size  distributions 
(Fig.  18)  at  Mace  Head,  on  March  10-11 ,  1988. 


SODIUM  CHLORIDE 


WAVELENGTH  (pi) 


REFRACTIVE  INDEX 


TOTAL  EXTINCTION  COEFFICIENT  ( 


0.55 

1. 548-1. 0xl0-7  i 

1.542-1. 0x1 0“7  i 

2.25x10  5 

0.6934 

2.27xlO~5 

1.06 

1.53  -1 .0x10” 7  i 

2.06xl0~5 

3.8 

1. 513-1. 0xl0-7  i 

5.39xl0~6 

10.6 

1.49  -l.OxlO-7  i 

1.32xlU“7 

0.55 

0.6934 

1.06 

3.8 

10.6 


AMMONIUM  SULPHATE 

1.52- l.OxlO-7  i 

1.52- l.OxlO  7  i 
l.51-2.4xl0~6  i 
1.56-.02  i 
1.99-.06  i 


2.26xlO“5 

2.26xl0-5 

2.05xl0-5 

7.19xl0-6 

1.025xl0-6 


The  calculations  were  made  for  the  visible  wavelength  of  0.55  um  together 
for  the  wavelengths  corresponding  to  the  laser  lines  0.6934,  1.06,  3.3  and 
10.6  pm.  Calculations  of  the  total  backscatter  coefficient  (m  ^  sr  for 
the  same  data  set  and  for  the  same  five  wavelengths  are  given  in  Table  5 
for  the  constituents  of  sodium  chloride  and  ammonium  sulphate. 

The  differential  extinction  coefficient  (per  m  per  micron)  is  plotted  in 
Fig.  19  as  a  function  of  particle  radius  based  on  the  averaged  particle 


Table  5 .  Backscatter  coefficient  (a  ^  )  for  visible  up  to  middle  1R 

wavelengths  based  on  Measured  particle  size  distributions 
(Fig,  lb)  at  Mace  Head ,  on  March  10-11 ,  1988 . 


WAVELENGTH  (ym) 

SODIUM  CHLORIDE 

REFRACTIVE  INDEX 

TOTAL  EXTINCTION  COEFFICIENT 
(m-^  sr-^) 

2.48xl0-6 

2.42xl0-6 

0.55 

1. 548-1. 0xl0~7  i 

0.6934 

1. 542-1. 0x10" 7  i 

1.06 

1.53  -l.OxlO-7  i 

1. 513-1 .0xl0-7  i 

1.49xl0-6 

1.025xl0-7 

3.8 

10.6 

1.49  -1.0xl0~7  i 

l.69xl0-8 

AMMONIUM  SULPHATE 

0.55 

1.52-l.OxlO-7  i 

l.SGxlO-6 

0.6934 

1.52-l.OxlO-7  i 

1.98xl0~6 

1.06 

1.51-2.4xl0-6  i 

1.20xl0-6 

3.8 

1.56-.02  i 

1.017xl0-7 

10.6 

1 .99-. 06  i 

5.31xl0-3 

size  distribution  data  of  four  cycles  for  'larch  10th,  1805  through  2120 
(shown  in  Fig.  14).  The  particles  are  assumed  to  be  comprised  of  ammonium 
sulphate.  The  calculations  are  for  wavelength  0.55  pm  (uppermost  curve) 
through  to  10.6  pc  (lowest  curve).  A  similar  plot  is  show  in  Fig. 20  for 
sodium  chloride  particles.  Figures  21  and  22  show  calculated  values  of 
the  differential  backscatter  coefficient  (per  .n  per  steradian  per  micron) 
as  a  function  of  particle  radius  for  the  same  data  set. 

Calculated  values  of  the  differential  extinction  coefficient  (per  m  per 
micron)  as  a  function  of  particle  radius  are  shown  in  ’-ig.  23  and  24  for 
ammonium  sulphate  and  sodium  chloride  particles  based  on  the  averaged 
values  of  particle  size  distribution  measured  over  five  cycles  for  the 
March  lOth-llth,  1988,  2230  through  0235  (see  Fig.  IS).  Figures  25  and  26 
give  the  differential  backscatter  coefficient  as  a  function  of  particle 
radius  for  the  same  data  set . 

Finally,  the  extinction  to  backscatter  ratio  calculated  for  each  of  the 
particle  probe  size  channels  at  wavelengths  0.55,  0.6934,  1.06,  3.8  and 
10.6  pm  are  plotted  in  Fig.  27  and  Fig.  28  for  constituents  ammonium 
sulphate  and  sodium  chloride.  One  can  see  that  the  extinction  to 
backscatter  ratio  approach  constant  values  for  the  visible  to  near  infrared 
wavelengths  but  have  resonant  values  at  the  larger  middle  IR  wavelengths. 
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Time:  0330-0911.  Wind  direction:  easterly 


icron  as  a  function  of  particle  radius,  (uncorrected) 


il  particle  concentration  per  c.c  per  in 
lead,  11  Dec.  1987.  Time:  0330-09 


cc  as  a  tunction 


Fig.  7.  Aeros 


osol  particle  concentration  per  cc  per  micron  as  a  function  of  particle  size 
e  Head  14  January  1988.  Time  1140-1145.  Wind  direction:  south  westerly. 
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Fig.  9.  Aerosol  particle  concentration  per  cc  as  a  function  of  particle  radius. 

Mace  Head,  15  January  1988.  Time:  0038-0043.  Wind  direction:  south  westerly 


Four  consecutive  data  cycles  of  particle  number  concentration  per  cc  as  a  function 
of  particle  radius. 

Mace  Head,  10  March  1988.  Time:  1805-2120.  Wind  direction:  westerly. 


Four  consecutive  data  cycles  of  particle  number  concentration  per  cc  per  micron 
as  a  function  of  particle  radius. 

Mace  Head,  10  March  1938.  Time:  1805-2120.  Wind  direction:  westerly. 
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Fig.  15.  Five  consecutive  data  cycles  of  particle  number  concentration  per  cc  as  a  function 
of  particle  radius. 

Mace  Head,  10-11  March  1988.  Time:  2230-0235.  Wind  direction:  west  south  westerly. 


Fig.  16.  Five  consecutive  data  cycles  of  particle  number  concentration  per  cc  per  micron  as 
a  function  of  particle  radius. 

Mace  Head,  10-11  March  1988.  Time:  2230-0235.  Wind  direction:  west  south  westerly. 
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rage  of  the  cyclic  data  shown  in  Fig.  16. 
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Differential  extinction  coefficient  per  m  per  micron  as  a  function  of 
particle  radius  and  wavelength  based  on  averaged  particle  size 

distribution  for  the  period  1ROS-2120,  March  10th,  1988.  Constituent:  Ammonium  sulphate. 
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Fig.  20.  Differential  extinction  coefficient  per  m  per  micron  as  a  function  of 

particle  radius  and  wavelength  based  on  averaged  particle  size  distribution 
for  the  period  1805-2120,  March  10th,  1988.  Constituent:  Sodium  chloride 
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Differential  backscatter  coefficient  per  m  per  sr  per  micron  as  a  function  of  particle 
radius  and  wavelength  based  on  averaged  particle  size  distribution  for  period 
1805-2120,  March  10th,  1988.  Constituent:  Ammonium  sulphate. 
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Fig.  22.  Differential  backscatter  coefficient  per  m  per  sr  per  micron  as  a  function  of  particle 
radius  and  wavelength  based  on  averaged  particle  size  distribution  for  period 
1805-2120,  March  10th,  1988.  Constituent:  Sodium  chloride. 


1.0- 


Di f f erential  extinction  coefficient  per  m  per  micron  as  a  function  of  particle 
radius  and  wavelength  based  on  averaged  particle  size  distribution  for  period 
2230-0235,  March  lOth-llth,  1988.  Constituent:  Ammonium  sulphate. 
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Fig.  25.  Differential  backscatter  coefficient  per  m  per  sr  per  micron  as  a  function  of  particle 
radius  and  wavelength,  based  on  averaged  particle  size  distribution  for  period 
2230-0235,  March  lOth-llth,  1988.  Constituent:  Ammonium  sulphate. 
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Fig.  26.  Differential  backscatter  coefficient  per  m  per  sr  per  micron  as  a  function  of  particle 
radius  and  wavelength,  based  on  averaged  particle  size  distribution  for  period 
2230-0235,  March  lOth-llth,  1988.  Constituent:  Sodium  chloride. 
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Fig.  28.  Extinction  to  backscatter  ratio  per  sr  as  a  function  of  particle  radius 
and  wavelength.  Constituent:  Sodium  chloride. 


